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Abstract 



Mh ' In the present work we undertake a re-examination of effective mass theory 

(EMT) for a semiconductor quantum dot. We take into account the fact that 



the effective mass (rrii) of the carrier inside the dot of radius R is distinct from 
the mass (m ) in the dielectric coating surrounding the dot. The electronic 



en 

structure of the quantum dot is determined in crucial ways by the mass dis- 
On 



continuity factor j3 = mi/m . In this connection we propose a novel quantum 



scale, a, which is a dimensionless parameter proportional to (3 2 VoR 2 , where Vq 



o 

represents the barrier due to dielectric coating. The scale a represents a mass 
modified "strength" of the potential. We show both by numerical calculations 
and asymptotic analysis that the charge density near the nanocrystallite sur- 



/\ • face, p(r = R), can be large and scales as 1/a. This fact suggests a significant 



role for the surface in an EMT based model. We also show that the upshift 
in the ground state energy is weaker than quadratic, unlike traditional EMT 
based calculations, and chart its dependence on the proposed scale a. Finally, 
we demonstrate that calculations based on our model compare favorably with 
valence band photoemission data and with more elaborate theoretical calcu- 
lations. 
PACS INDEX: 71.24.+q, 73.20.Dx, 78.66.-w, 61.46,+w 
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I. INTRODUCTION 

Semiconductor nanocrystallites, more popularly known as quantum dots (QDs), have 
been extensively studied in the past decade and a half. The system is interesting from 
the point of view of basic physics, with the carriers being confined to an essentially "zero" 
dimensional structure. The efficient luminescence observed in some of these crystallites 
makes them promising candidates for opto-electronic devices. Further, the inexorable drive 
towards device minituarization makes them technologically significant. 

The physical dimensions of a typical QD are in the range 1-10 nm. It is coated by a 
dielectric which may be a polymer, glass, oxide, etc. depending on the method of prepara- 
tion. In a simplified effective mass theory (EMT) approach, the QD is taken to be spherical 
with radius R and the dielectric presents a finite barrier Vq 6 [1 - 10 eV] to the carriers 
(electrons/holes). This is shown in Fig. |l|. The carrier confinement leads to an enhancement 
of the "band" gap which, in the simple quantum confinement model (QCM), scales inversely 
with the size (1/R 2 ). However, experimental observations and more elaborate theoretical 
calculations suggest that the gap dependence is infra-quadratic (1/-R 7 , 1 < 7 < 2). Further, 
several workers have questioned QCM and suggested that it downplays the significance of 
surface-related effects. 

In the present work we demonstrate that EMT can be reconciled with the observed 
infra-quadratic shift of the "band" edges. More importantly, an examination of the charge 
density in our EMT based model reveals the significant role of the surface in a direct fashion. 
To accomplish this, we re-examine the EMT taking into account the fact that the effective 
mass of the carrier inside the dot (rrii) is distinct from the mass outside (m ). In other 
words, there is a mass discontinuity across the potential barrier depicted in Fig. [l|. The 
ratio (3 (= rrii/m ) can significantly affect the electronic properties of a QD. It is well known 
that the classical scale representing the surface to volume ratio (S/V) plays an important 
role in determining both classical and quantum properties. We propose a novel quantum 
scale a which is a dimensionless parameter proportional to j3 2 V R 2 (Eq. (|15|)). This scale 



represents a mass modified "strength" of the potential (MSP). The correction to the ground 
state energy E 1 and the charge density at the interface p{r = R) scale in an appealingly 
simple way with the proposed scale a. 

In Sec. we describe the basic model in brief. The carrier (electron/hole) is confined to 
a spherical QD of radius R with a finite barrier of height Vq (see Fig. [I]). As suggested by 
BrusEl, the BenDaniel-Duke boundary conditionsa are appropriate for a system with a mass 
discontinuity across the barrier. This aspect has been ignored in a number of theoretical 
calculations in the past (see Sec. [V|). As will be shown subsequently, the physics of the 
problem is determined in crucial ways by the mass discontinuity factor /3(= mi/m ). A 
pedagogical point, we highlight is that the BenDaniel-Duke boundary conditions must be 
applied to the full wavefunction and not to the partial wavefunction. The error involved in 
employing the partial wavefunction can be substantial and is presented in Table |. 



The numerical results are presented in Sec. |T|. Previous works in the field have devoted 
considerable attention to the upshift of the ground state energy E\. We focus instead on the 
charge density p(r). BrusEl had pointed out that the charge density near the nanocrystallite 
surface can be large if the carrier effective mass rrii is small. We show that the charge density 
at the nanocrystallite boundary p(r = R) is indeed large if (3 (= mi/m ) is small (Fig. ^. 
A more complete picture is provided by our proposed MSP scale a (oc (3 2 VqR 2 ). In Fig. |] 
we demonstrate that p(R) scales inversely with a . For a quantum dot, S/V ratio is large. 
Thus the fraction of unsaturated bonds on the surface increases and/or the formation of 
molecular complexes on the surface is favoured. The large charge density on the surface 
suggests a novel non-classical reason for the importance of the surface. We also show that 
the ground state energy E\ exhibits a size dependence which is weaker than the quadratic 
(1/i? 2 ) dependence one has come to expect from EMT (Fig. |5]). 

A detailed asymptotic analysis is presented in Sec. [TV] to explain the numerical results of 
Sec. [T IT] and to present them in terms of the proposed quantum scale a. We show that p(R) 
scales as ~ 1/er. We introduce the notion of the peneteration depth S in this connection. 
We also explain the infra-quadratic dependence of the ground state energy on the size R. In 



particular, we demonstrate that the correction to the quadratic term scales as ~ —\j\fo. 

Section [V] constitutes the discussion. We show that our model compares favourably 
with the experimental valence-band photoemission data of Colvin et. a/.S on CdS (Fig. 0) 
and with more elaborate tight-binding calculations^ (Fig. |7|). We summarize our work and 
suggest directions for future research. 

II. BASIC THEORY 

We consider a spherical semiconductor nanocrystallite of diameter d = 2R (See Fig. HJ). 
An additional electron in the "conduction band" of such a crystallite is described in effective 
mass theory by the Hamiltonian: 

H = ~^V • f-i^- • v) + V(f) (1) 

2 \m*(r) J 

For a position dependent mass the appropriate hermitian kinetic energy operator is given 
by the first term on the right hand side of Eq. (0)0. The electron effective mass inside the 
nanocrystallite {mi) is different from the effective mass in the dielectric coating (m ). A 
useful parameter is the ratio (f3) of the effective masses. 

P = — 2 

m 

For the experimentally relevant cases explored in the last two decades, (3 is less than unity. 
For extreme cases, such as GaAs or InAs quantum dots, /3 equals 0.067 and 0.02, respectively, 
where we have taken m = m e , the free electron mass. An important aspect of the present 
study is to highlight the significance of this parameter (3. The potential V(r) in Eq. ([[]) is 
given by 

{0 r < 

(3) 
V r> R 

Here V is a large positive potential and represents the dielectric coating surrounding the 
nanocrystallite (See Fig. |l|). Typically a nanocrystallite is surrounded by dielectrics such 



as glasses, polymers, organic solvents or oxides and hydridesH. Thus, the electron is in a 
spherically symmetric well. The wavefunction for this spherically symmetric problem can 
be written in the form: 

^ n i m {r) = Rni{r)Y lm {e,4>) (4) 

where the symbols have their usual meanings. It is customary to write the radial wavefunc- 
tion R n i as: 

u n i{r) = rR n i (5) 

The equation for the partial wavefunction u n i(r) is : 

d 2 u nl 1(1 + 1) 2rrnE 

-&T ~ —^uni + —u nl = (6) 

where rrii is the mass of the electron inside the potential well. The electron effective mass 
m* assumes different values mj inside and m outside the well. 

For the I = case, the standard form for the radial wavefunction is obtained on solving 
Eq. (||). This is : 

*»M=A*g£> (7) 



n 



for < r < R and 



R no(r) = B— (8) 



for r > R, where 



C = f^ 0) 



and, 



V 7/ 

and A and B are normalization constants. 



out _ , / 2m (K - E n ) (1Q) 



Standard text books on quantum mechanics state the condition of continuity of the 
derivative of the wavefunction inside and outside the well as : 



dR 



ni 



dr 



dR 



nl 



t>- dr 



(11) 



However, as has been discussed in the semiconductor literature, this condition must be 
replaced by the BenDaniel-Duke condition in case the effective masses are different across 
the interface. The condition now reads 



1 dR 



nl 



raj dr 



1 dR, 



ni 



r ^ R - rrio dr 



(12) 



A further point needs to be made at this stage. Many textbooks, in the process of the 
solution of the Schrodinger equation, define the partial wavefunction u n i (Eq. fl5|)) and then 
further state that the continuity conditions at the interface may be imposed on u n i(r)uu. 
They claim that this would yield the same eigenvalue conditions as imposing these conditions 
on R n i(r) would. Quoting Schifip, " The energy levels are obtained by ...(this is equivalent 
to making 1/R dR/dr continuous there)". It is shown that this is manifestly incorrect when 
we use the modified boundary condition. It is further shown that in the limit the effective 
masses become equal across the interface, the two treatments indeed become equivalent for 
this special case. Imposition of the modified boundary condition on the full wave function 
leads, after some algebraic manipulations, to 



urn* , l ~P 



kZcot{kZR)=-pkZ* + - 1 f. (13) 

whereas, the use of the same condition on the partial wavefunction u n i(r) leads to 

ccot(Ci?) = -/%r (14) 

The two conditions are evidently quite different. However, they agree when f3 is unity, recov- 
ering thereby the elementary case. There are three parameters, j3, Vq and R, in the above 
transcendental equation. As can be readily seen from the two equations, the eigenvalues 
calculated using Eq. fll3|) would deviate more and more from those calculated using Eq. fll4l) 
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for smaller (3 and smaller R. The energy eigenvalues obtained for a particular set of param- 
eters is presented in Table Q. For /3=0.1, Vo=1.0 and i?=20.0 A°, the disagreement is even 
worse, being as large as 57%. In the sections that follow, we shall employ Eq. ([13]). 

III. RESULTS 

In this section, we present the results of our numerical calculations. Further, this section 
and the next contains the interpretation of our results. In what follows, we focus our 
attention on the spherically symmetric (I = 0) case. 

A quantity of central interest to us is the radial charge density p{r) 

p{r) = r 2 R 2 n0 (r) 

Where the multiplicative constant, e, the carrier charge, has been ignored. We examine 
the dependence of the equilibrium charge density on the ratio of the effective masses (3 
(= rrii/m ), the barrier height Vq and the radius of the nanocrystallite R. In Fig. |2|, we 
have plotted the charge density for different values of (3 and Vq — 1.5 eV. The radius of the 
crystallite is R = 50 A°. As can be readily seen, the value of the equilibrium charge density 
at the crystallite boundary, p(r = R), falls rapidly with increasing (3. The peak in the charge 
density occurs at r G [R/2,R]\ in other words, close to the surface of the crystallite. For 
small values of /3, the peak is at r ~ R. Examples of small (3 are InP, InAs, GaAs, etc. Even 
if the carrier is photogenerated in the interior of the crystallite, it rapidly redistributes itself 
such that the equilibrium charge density is large near the boundary. In nanocrystallites with 
large (3{~ 1) the peak in the charge density shifts towards R/2. In other words r = R/2 is 
a "fixed point" for the charge density. Figure ^ shows this more systematically. We have 
plotted the charge density at the boundary as a function of (3, and taking Vo = 5 eV. There 
is a clear decrease in the charge density as the value of j3 increases. In the next section, 
we shall show this more rigorously in the limit of large Vq. An appropriate perspective on 
p{r = R) is obtained by defining a mass modified " strength" of the potential (MSP), a. 



a = (fiKoRy (15) 



Kq 



'2m Vo 



I h 2 

Conventional textbooksBEI define the strength of potential as VqR 2 . As will become apparent 
in the next section, a more appropriate definition for the case under study is a = (/3kqR) 2 . 
Figure |3] depicts that p(r) falls with o. This decrease is observed to be nearly linear. We 
shall have occasion to examine this universal parameter parameter a in the next section. 
In Fig. H we study the dependence of the ground state energy E\ (n — 1, I — 0) of the 
carrier on the crystallite size R. We have taken Vq = 5 eV for this figure. We observe that 
this dependence is infraquadratic, i.e. weaker than 1/R 2 . We note that this infraquadratic 
behaviour has been reported both experimentally and on the basis of tight-binding (TB) 
calculationstyty. In the present case infraquadratic behaviour can be traced to the finite 
nature of the barrier and the change of the effective mass across the barrier. Inspection of 
Fig. |5] reveals that with decreasing value of /3, the behaviour increasingly departs from the 
quadratic case. This is demonstrated rigorously in the next section on asymptotic analysis 
(Eq. (|23[)), where we see that 

C\ and C2 being constants. An attempt to coerce the above expression into the form 

E, - % (17) 

will yield 1 < 7 < 2. We have presented the salient results of our model in this sec- 
tion. Additional results related to experiments and earlier theoretical calculation on specific 
semiconductor nanocrystallites will be presented in the last section (Sec. 0). The following 
section on asymptotic analysis will attempt to explain the numerical results obtained in this 
section. 



IV. ASYMPTOTIC ANALYSIS 



We now present an asymptotic analysis of Eq. (O) in an attempt to explain the results 
of our calculations in the previous section. For an infinite potential well, the eigenvalue 
condition reduces to : 



kTR = nn 



(18) 



where k™ is defined in Eq. ([]). We consider the case of the ground state i.e. n = 1. For a 
well which is "sufficiently" deep but finite, we approximate the above expression by : 



k\ n R = 7i - e 



(19) 



where e is a small number. Using Eq. 



in Eq. ([13|) , we obtain 
1 



TV 



where 



13(1 + k R) 



l 2m V 



hf 



(20) 
(21) 



where we have assumed that V ^> Ei, the ground state energy. The condition on the 
smallness of e is now apparent. If either of these parameters (3, R, V is large then e would 
be small. It is clear that e is inversely dependent on the mass modified "strength" of the 
potential defined earlier in Eq. (|T5| ) 

e c * (22) 

V a 

when yfo ^> j3. We are now in a position to obtain an approximate expression for the 
ground state energy E\. Using Eqs. flT9| ) and ( |20|) and the definition of k™ from Eq. ([|), 
one obtains, 



2*2 



Ex 



TX l % 



2 mi R 2 

n 2 h 2 
2m,:R 2 



1- 



f3n R 



+ . 



1- 



+ 



(23) 
(24) 



Thus one can see that the ground state energy has the size dependence depicted in Fig. |5|, 
namely, 

- 15 < 26 > 

with effective exponent 7 < 2, as mentioned earlier. The EMT literature on nanocrystallite 
semiconductors commonly quote 7 = 2. However, absorption and luminescence experiments 
as well as the tight-binding calculationstiM yield 7 < 2 . We have thus formally demon- 
strated how even within the EMT approach, an infraquadratic exponent (7 < 2) is obtained. 
The calculation depicted in Fig. |] attests to this. 

Garrettll3 has proposed a penetration depth for the one dimensional finite well problem 
with (3=1. We define an analogous penetration depth (5) for the carrier in the three 
dimensional potential well and for the general case /3 ^ 1. 

/3k 
" (27) 



a 



Next we shall examine the charge density at the surface. An appealing, intuitive way to 
understand the nature of p{r = R) is to use the concept of the penetration depth. It is 
easy to see from Eq. ( p7[ ) that the penetration depth will be large in case f3 is small. This is 
already a confirmation of the result that the charge density is large at the boundary for small 
/3 (Figs. |2| and ^). Thus one would expect that when the carrier mass in the semiconductor 
is small (InSb, GaAs, CdS) then the charge density at the interface will be large. The 
penetration depth 8 indicates the extent of the wavefunction penetration into the forbidden 
region. 

We can present a more rigorous analysis of p(r = R) by examining the wavefunction. 
The normalized wavefunction inside the potential well is : 

10 



where 



.4 



k[ n 



R ( _ sin(2Jfej n .R)\ 1 - cos(2k[ n R) 
~2 I 2k\ n R J + 4fcp 



(29) 



Using Eq. C [T9|) in Eq. (]2q) , we obtain after some algebraic manipulations: 



R 10 (r) « A 1 + - 



7T 



sin(7rr/i?) e /7rr 
vrr/i? ~^ COS l^ 



(30) 



whence we get for the charge density p, 



p(e,r) = r 2 R 2 w (r)=A 2 (l 



7V 



R\ 2 9 (ixr\ fre\ 2 9 / ;r/ 
) sin 2 ( — + — cos 2 



TV J 



RJ 



tv / 



R 



(31) 



The charge density at the interface is: 

p(R) = r 2 Rl (r) 



r=R 



A 2 (l + -)(-) R 2 

71/ \7T 



(32) 



Recall that e 2 ~ 1/a (= 1/(/3k R) 2 ). It is therefore clear that p(r = R) would fall with 
increasing j3. This prediction is in agreement with Fig. [| where p(R) is l//? 185 . 

Within the range of validity of approximation in Eq. (|19D , e is much less than unity. 
Therefore, Eq. (|32|) may be written as 

,R 2 



p(R) w A 2 



a 



(33) 



One can thus understand the inverse dependence of p(R) on a in Fig. f|. It is interesting to 
note that the charge density can be related to the penetration depth defined in Eq. (|27|). 



2x2 



p(R) w A 2 S 

One may also obtain the peak position, r pea k, in the charge density by a simple differentiation 
of ri?io(r), where -Rio(^) is given by Eq. (|30|) . This yields 



(nr peak \ (Trr peak \ 1-e/vr 

l-R-J tan l-R-J = -^r 



(34) 



It can be seen easily, either by plotting Eq. (|3~TD or by solving Eq. ( |3^ ) numerically, that the 
peak in the charge density p is located at values of R/2 < r < R and shifts to r = R/2 as 
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the value of e is decreased. Decrease in e is brought about by an increase in a. Equation fl34|) 
indicates that r peak e [R/2,R]. When (3 is small, r peak > R/2 and this is also graphically 
borne out by Fig. 0. Since the peak is shifted towards the crystallite boundary, the charge 
density at the boundary becomes substantial. When the MSP scale a is large, Eq. ( |34| ) 
implies that r peak — >■ R/2. In other words there is a fixed-point or "attractor" at r = R/2. 
The peak in the wavefunction does not shift to a value below R/2, no matter how large a 



is. 



V. DISCUSSION 

The objective of the present work is to explore EMT in the context of the variation of 
the carrier effective mass across semiconductor - dielectric interface in a QD. We consider a 
simple Hamiltonian described by Eq. ([I]) in order to illustrate the essential physics. It would 
be worthwhile to compare our results with other theoretical calculations and experiments. 
Below we present two such comparisons. 

We compare our results with the valence band photoemission experiment carried out 
by Colvin et. al. on CdS quantum dotsS. Their work represents the first non-optical 
observation of the electronic structure of semiconductor nanostructures. The experiment 
measures the energy shift of the valence band maximum as the cluster size decreases, taking 
the largest cluster of R = 35 A" as a reference. For this caseQ, j3 = 0.53 and Fig. |] depicts 
the dependence of the valence band shift on the cluster size as calculated by us. We fit 
our data to C/R a . We get an exponent a = 1.26, which is in good agreement with the 
experimental observation. Colvin et. al. attribute this shift to two factors : (i) the kinetic 
energy enhancement due to the quantum confinement of the hole ; (ii) the polarization of 
the crystallite and the loss of dielectric solvation energy™!^. As is clear from Fig. |^ a finite 
barrier quantum confinement model with appropriate effective masses provides a reasonable 
explanation for the experimental behaviour. We also note in passing that Takagahara has 
argued that the polarization term can in general be neglectedE-a. 
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Lippens and LannooQ have theoretically studied the same system. Our model can be 
compared with their results for the valence band shift. Fig. ^ depicts such a comparison. 
We have taken the barrier height Vo = 2 eV. A majority of the literature claim that TB 
calculations are in better agreement with experiments. The EMT with quadratic dependence 
on size (~ 1/-R 2 ) represents an extreme case. As can be seen from Fig. [!] our EMT based 
model is in close agreement with the experimental observations. On the other hand, the 
standard EMT calculation of Lippens and Lannoo is at variance with their TB calculation. 

The majority of the theoretical work for the absorption and luminescence in quantum 
dots takes excitonic effects into account. In other words, the Hamiltonian involves both a 
hole and an electron. We propose to extend our simple EMT model to the excitonic case as 
well as to study excited states in the future. We also plan to study the capacitance of QDs 
by analyzing multi-electron effects in our model. 

We emphasize that the EMT calculation should be carried out with the correct boundary 
condition as presented in Eq. flT2|). This fact has been pointed out by Brusa. The use of this 
boundary condition leads to an eigenvalue condition (Eq. flT3])) which is different from the 
one normally encountered in textbooksQtl (Eq. (|TJ]) with (3=1). A further point to note is 
that the boundary condition must be imposed on full wavefunction R n i(r). Imposing it on 
the partial wavefunction u n i(r) would lead to incorrect results. 

EMT based calculations have normally been carried out using an infinite barriered. In this 
case, the charge density at the nanocrystallite boundary is zero. It is therefore not surprising 
that the significance of the surface is downplayed. Further, workers who recognize the finite 
nature of the barrier have ignored the mass discontinuity across the barrierll3il3 and the 
reason sometimes stated is that most of the population density of the carrier is still confined 
well within QDli3. The LDA based calculations on dot capacitance and shell filling effects also 
ignore the mass discontinuity^. Porous silicon, which is a disordered agglomeration of silicon 
nanocrystallites, presents a case-study for the debate between the quantum confinement 
and the surface state schools of thought. This debate has been reviewed by a number of 
workerstirEil. In this work we have shown that the charge density can be large at the surface 
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even in an EMT based framework. This is more so for semiconducting materials with small 
carrier effective masses and dielectric coatings representing a small potential barrier. Hence, 
the contention that surface related effects are significant appears valid. 

Another demerit of the infinite barrier model is its prediction of an inverse quadratic 
(1/R 2 ) shift of the band gap. We have already stated that this does not agree with exper- 
iments and with more elaborate theoretical calculationsliyUc3'c3. We have shown through 
our work that a calculation for finite barrier model leads to a good agreement with the 
experimental infraquadratic dependence. 

The highlights of the work are : 

1. Ground state energy within EMT scales infraquadratically with the crystallite size 
(Eq. (H) - ©). 



2. The parameter of relevance is our proposed mass modified "strength" of the potential, 
the MSP scale a = (/3k R) 2 . Our proposed parameter a presents a more complete 
picture of the physical situation thau (, aloue doesl 

3. The charge density at the boundary p(R) can be large (Figs (Q) - (ffl)) and scales in- 
versely with a. This is unlike what one normally expects from EMT. In the literature 
of semiconductor nanocrystallites there has been a longstanding debate between the- 
ories for quantum confinement versus surface statescJlHJ. The present work provides 
a fresh perspective on the importance of surface related phenomena. 
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TABLES 



TABLE I. Comparison of eigenvalues obtained by applying the correct boundary condition on 



the full wavefunction (Eq. (|13|)) and on the partial wavefunction (Eq. (14)). The values of the 
parameters employed are /3 = 0.1, Vq = 2.5 eV and R = 40 A°. 



Eigenvalues with 
full wavefunction 


Eigenvalues with 
partial wavefunction 


Percentage difference 


0.13 


0.15 


12.04 


0.61 


0.65 


3.21 


1.52 


1.56 


2.53 
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FIGURES 
FIG. 1. The left side of the figure depicts an idealized spherical quantum dot (QD) surrounded 

by a dielectric coating. This is modeled by a potential well of height Vq. The carrier effective mass 

is mi inside the well and m in the dielectric coating outside. 

FIG. 2. A typical normalized charge density p(r) inside a crystallite of size R = 50 A. Note 
that for small fi = (mi/m ), p(r) is large at the crystallite surface (r = R). As (3 increases, the 
position of the peak r pea f. in p shifts to R/2. 

FIG. 3. The charge density at the crystallite boundary, p{r = R). The charge density p(R) falls 
with increasing (5. We fit our data for p{R) to C\j ' (3 a , with a = 1.85. This is in good agreement 
with the prediction of Eq. (|32j). C\ is a constant. 

FIG. 4. The dependence of the charge density p(R) at the crystallite boundary on the proposed 
quantum scale a (= ((3kqR) 2 ). The scale a represents a mass modified "strength" of the potential 
(see text). The data for p{R) fit well to C^/o"" with a = 0.98, and in agreement with the prediction 
of Eq. d33|). C2 is a constant. The dotted line is a fit to the calculated data (diamonds). 

FIG. 5. The ground state energy E\ versus the crystallite size R. The calculated data are for 
four different values of (3 and are fitted to C/R 1 , where C is a constant. The values of 7 for the 
different values of (3 = 2.0, 0.4, 0.1 and 0.01 are respectively 1.95, 1.76, 1.39 and 1.03. We can see 
that the energy shift is increasingly infraquadratic for smaller values of /?. 

FIG. 6. The data taken from the valence band photoemission spectra for CdS crystallite due 
to Colvin et. al. (Ref. 3) is indicated by diamond symbols. Our calculation is indicated by the 
dotted line. A fit to C/R a leads to a = 1.26 for both cases. Our calculation gives a much improved 
result as compared to traditional EMT with a = 2. See Sec. for further discussion. 



FIG. 7. Valence band shift with size of the CdS crystallite. Our results (solid line) are very close 
to the tight-binding (TB) calculations (diamonds) of Lippens and Lannoo (Ref. 4). In contrast, 
there is less agreement between the EMT (dashed line) and TB (diamond symbols) calculations of 
Ref. 4. 
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